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Abstract 

We report on first-principles total-energy and phonon calculations that clarify structural stability 
and electronic properties of freestanding bilayer silicene. By extensive structural exploration, we 
reach all the stable structures reported before and find four new dynamically stable structures, 
including the structure with the largest cohesive energy. We find that atomic protrusion from 
the layer is the principal relaxation pattern which stabilizes bilayer silicene and determines the 
lateral periodicity. The hybrid-functional calculation shows that the most stable bilayer silicene is 
a semiconductor with the energy gap of 1.3 eV. 
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Bilayer graphene provides new aspects of graphene physics such as the band gap 
opening-"- and moire-pattern-induced electron localization.-!^ The interlayer interaction 
is obviously weaker than the intralayer interaction but decisive to modulate the electron 
states due to its symmetry breaking. Similar intriguing behavior with the extension related 
to the spin-degrees of freedom is expectect^"- for layered Si (silicene), which has been grown 
experimentally with a form of monolayer—"— and of a few layers!^"— mainly on Ag sub¬ 
strates. One of the most important characteristics which discriminate silicene from graphene 
is the buckling of two sublattices caused by the preference of Si for sp^ hybridization. In 
fact, first-principles calculations within the local density approximation^^i^i have clarified 
that a planar Si monolayer is unstable to the buckling, and that the resultant freestanding 
monolayer silicene with the buckling of 0.44 A have the Dirac cone at the Fermi level, 
Ep. This buckling brings about a complex but rich variation in structure and electronic 
properties in bilayer silicene, that we address in this Rapid Communication. 

Several theoretical investigations have been performed on freestanding bilayer silicene.—"— 
A fiat bilayer structure with the perfect overlapping stacking (AA stacking) has been pre¬ 
dicted based on a molecular dynamics simulation.— A corrugated 2x2 (rectangular su¬ 
percell) reconstructed structure different from the well known yr-bonded chain structure of 
the Si(lll) surface^^ has been also foundRecently, freestanding bilayer silicene has been 
found to have several local minima in total energy as a function of the lateral lattice parame¬ 
ter, suggesting its wealth of the structural diversity.— A possible structural phase transition 
under the lateral strain has been also discussed.— To validate such theoretical predictions, 
however, a systematic exploration of stable structures with their thermal excitation spectra 
(phonons), which is lacking in the past, is imperative. Further, clarification of correlation 
between the structural diversity and its role in electronic properties is highly demanded. 

We here perform systematic first-principles total-energy and phonon calculations for free¬ 
standing bilayer silicene with various lateral periodicities and atomic densities. By extensive 
geometry optimization followed by the phonon calculations, we unequivocally identify ten 
dynamically stable structures with distinct atomic configurations, symmetries, and period¬ 
icities. We clarify that the ten structures include all the six structures reported in the past. 
Other four structures are newly found and more stable than the previously reported ones. 
We find that the more stable structures have a single prominent structural characteristics, 
i.e. the protrusion of Si atom. We also find that there is an energetically optimum lateral 
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periodicity for the protruded structure, i.e., ^/3 x ^/3 or 2 x 2, depending on the stacking 
of the two Si layers. We further clarify that the stable freestanding bilayer silicene is a 
semiconductor in which the energy gap is sensitive to the detailed protruded structure. 

We use the pseudopotential-planewave method2i>^ based on the density-functional 
theory.—!^ We adopt the exchange-correlation functionals developed by Perdew, Burke, and 
Erenzerhof (PBE).— We crosscheck the electronic properties by using the Heyd-Scuseria- 
Ernzerhof (HSE) hybrid functional.— >2^ Phonon dispersion relations are calculated based 
on the density functional perturbation theory.—"— Computational details are described in 
Supplemental Material.— 

In our calculations, we have considered 1 x 1, 2 x 2, \/3 x y/3 and y/T x ^/7 lateral period¬ 
icities with respect to the monolayer periodicity and performed total-energy optimization. 
Ereestanding monolayer silicene is intrinsically buckled so that the stacking of two Si mono- 
layers is of rich variety: the buckling of two layers may be in-phase, out-of-phase, or planar 
both in the AA stacking and in the AB stacking (Bernal stacking named in graphite); fur¬ 
ther the stacking is not restricted to AA or AB but could be orthorhombic (OR) in general. 
We have then noticed that there are 16 distinct stacking configurations.— Eor each stacking 
configuration, we have adopted a certain lateral lattice parameter and fully optimized atomic 
geometries, and then repeated the calculations with varying the lattice parameter. We have 
also examined a possibility of the unit-cell distortion. By such extensive exploration, we have 
obtained 24 distinct total-energy minimized geometries for the 1x1, the 2x2 and s/S x y/3 
periodicities. We have then calculated the phonon spectra for thus obtained total-energy 
minimized structures. Surprisingly, the 14 of the 24 total-energy minimi zed structures is 
clarified to be unstable with imaginary phonon frequencies, although their interlayer bind¬ 
ing energies are positive. We have then eventually reached four 1x1, four 2x2 and two 
\/3 X \/3 dynamically stable structures. Details of our structural exploration are described 
in Supplemental Material.^ 

The five dynamically stable structures obtained in the present study shown in Eig. S4 
(Supplemental Material^) and Eig. [11(a), which are named AA-lxl^ AB-lxl, slide-lxl, 
OR-1 xl and slide-2 X 2, are presumably identical to the structures which have been reported 
previously.—"— The AA-lxl structure is similar to the AA-stacking bilayer graphene, but 
the interlayer distance is significantly small (2.411 A), resulting in the fourfold coordination 
of all the Si atoms. The AB-lxl structure is similar to an atomic slab in the diamond 


3 


structure, inferring the existence of threefold coordinated Si. The slide-lxl structnre also 
shows threefold and fonrfold coordinated atoms, althongh the interlayer distance becomes 
shorter compared to the AB-lxl. The OR-lxl structure is obtained by relaxing the lateral 
and in-plane bond angles of AA-lxl structure. The AA-lxl, AB-lxl, slide-lxl, OR-lxl, 
and slide-2X2 structures are presumably identical to the 1AA-, the 1AB-, the slide-2AA- 
strnctures in ref. l27|, the Si-Cmme strncture in ref. l29|, and the phase II strncture obtained in 


ref. 


28l. respectively. We have also obtained the structure proposed in ref. 


2a, here denoted 


by rect-OR-2x2 [Fig. S4(e)], by nsing a rectangular 2x2 supercell. The shde-2x2 and 
the rect-OR-2x2 structures exhibit peculiar structural characteristics: two Si atoms in each 
unit cell protrudes prominently, causing the local tetrahedral geometry. It is noteworthy 
that this protrusion renders those two structures lower in energy than other 4 structures 
(see below). 

In the 2x2 periodicity, we have found two new structures, named AA-2x2 [Fig. He)] 
and hex-OR-2x2 [Fig. H®)]- The latter is the lowest in energy among the structures ever 
reported (see below). The AA-2x2 structure shows common characteristics to the slide- 
2x2 strncture [Fig. H^-) ]; but the details are different. The AA-2x2 structnre possesses 
almost complete AA stacking whereas the slight dislodgment is observed in the shde-2x2. 
In the AA-2x2, two silicon atoms in each unit cell are highly protruded [red spheres in 
Fig. He)] whereas another two atoms are moderately protrnded (green spheres). In the 
hex-OR-2x2 strncture, on the other hand, red and green protrnded Si atoms in Fig. H©) 
form dimers. The hex-OR-2x2 structure is similar to the rect-OR-2x2 structure although 
the supercell is not rectangnlar but hexagonal. The fonr 2x2 strnctures, inclnding newly 
found two structures, are dynamically stable, as is evidenced from their calculated phonon 
spectra shown in Figs. Ha), He) and He), and Fig. S4(j). 

We have further explored dynamically stable structures with different lateral periodicity. 
Foensing the dynamically stable stacking configurations in the 2x2 periodicity shown in 
Figs. (D (a), (c) and (e), we have examined the -^/S x \/3 lateral periodicity with those 
stackings. We have found that the AA- and OR-\/3 x \/3 structures [Figs. H^) and H^)] 
are also dynamically stable whereas the slide--s/S x structnre [Fig. Hb)] is unstable, as is 
evidenced in onr phonon dispersion relations in Figs.Hb),Hb) andH^)- To our knowledge, 
the A A- and the OR-v^ x -^/S structures found here have not been reported in the past. 

The calculated total energies of freestanding bilayer silicene with various stackings and 
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FIG. 1. (Color online) Total-energy minimized structures of x and 2x2 freestanding 
bilayer silicene: (a) slide-2x 2, (b) slide-x (c) 74^4-2x2, (d) AA-^/3 x (e) liex-OR-2x2, 
and (f) OR-v^ X \/3 structures. Unit cells are represented by silver lines. Blue, red, and green 
spheres represent fourfold coordinated, highly-protruded, and moderately-protruded silicon atoms, 
respectively. The angle between two lattice vectors are slightly changed from 120° in the slide- 
2x2, slide-\/3 x \/3, hex-OR-2x2, and OR-\/3 x \/3 cases. The XCrySDen program is used for 
visualization of the atomic structures;^ 


periodicities are listed in TableHl The hex-OR-2x2 structure newly found here has the lowest 
energy, and the slide-2 x 2 follows with the total-energy increase of 6 meV/atom. The OR- 
\/3 x \/3 and rect-OR-2x2 structures are also close in energy (-1-9 and -|-10 meV/atom from 
the lowest energy, respectively). The AA-\/2 x \/3 and AA-2x2 are in the second lowest 
energy group (-1-31 and -|-41 meV/atom, respectively). The OR-, AA-, slide-, AB-lxl 
structures have relatively high energy (-1-60, -|-78, -1-112 and -|-131 meV/atom, respectively) 
compared to the \/3 x \/3 and 2x2 structures. We have also calculated the total energy of 
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FIG. 2. Phonon dispersion relations of (a) slide-2x2, (b) slide-x (c) AA-2x2^ (d) AA- 
X (e) hex-OR-2x2, and (f) OR-\/3 X structures (in the same order as Fig. [1]). Phonon 

dispersions crossing several M points, M(0.5,0), M'(0,0.5) and M"(—0.5,0.5) in units of recip¬ 
rocal lattice vector, which are inequivalent to each other in the low symmetry cases are plotted. 
The negative values in (b) mean imaginary phonon frequencies, indicating that the structure is 
dynamically unstable. Imaginary frequency branches (acoustic vibrational modes along the direc¬ 
tion perpendicular to the layers) appear around the F point in other cases due to the interpolation 
errors. We have confirmed that these imaginary frequencies become real when we directly compute 
the phonon frequency at specific q points, with sufficient vacuum space (16 A) of the supercell 
and cutoff energy (60 Ry). 


another structure named honeycomb dumbbell silicene,— constructed from a 3 x 3 monolayer 
silicene. This dumbbell structure is outside our configurational phase space for the structural 
search, although it has higher total energy than the value of the hex-OR-2x2 structure by 
206 meV/atom (VS x a/3 periodicity) or 60 meV/atom (2x2 periodicity). 

Our systematic total-energy and phonon calculations have clarified rich structural variety 
of bilayer silicene. At the same time, it has become clear that stable bilayer silicene is accom- 
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TABLE 1. Calculated total energies (in meV/atom) of various stacking geometries {AA, slide, OR, 
and AB) and lateral periodicities (1x1, \/3 x \/3, 2x2, and V7 x \/7) of freestanding bilayer 
silicene. The total energies are compared to that of freestanding monolayer silicene. We do not 
find a distinct AB-stacking structure for the \/3 x \/3 and 2x2 periodicities. 


Stacking 



Periodicity 


1x1 

V3 X V3 

2x2 

vTx x/7 

AA 

-183 

-230 

-220 

-217 

Slide 

-149 

-234 

-255 

-223 

OR 

-201 

-252 

—261 (hex) 

—251 (rect) 

-219 

AB 

-130 

- 

- 

- 


panied with a characteristic relaxation, i.e., the atomic protrusion. The ^/3 x y/3 structures 
and their structural analogues with the 2x2 periodicity have common stacking configura¬ 
tions and protruded structures but different total energies. This implies the existence of an 
optimum periodicity for the protruded structure for each stacking configuration. In order to 
clarify this point, we have performed total energy calculations for the \/7 x \/7 structures 
having the characteristic protruded pattern for the three stacking geometries (AA, slide and 
OR). The total energies of the fully optimized y/l x \/7 structures are shown in Table HI 
They are all higher than the corresponding total energies of the 2 x 2-periodicity cases. The 
total energies of the slide and OR stacking geometries become minimum with the 2x2 peri¬ 
odicity, whereas the total energy of the AA-stacking is minimum at the y/S x \/3 periodicity. 
We have now unveiled that the atomic protrusion is the principal relaxation pattern which 
stabilizes the bilayer silicene and induces a particular lateral periodicity. 

Figure |S] shows the top view of the OR-\/7 x \/7 structure obtained here. Red and 
green spheres represent dimerized protrusion which is common to other periodicities. We 
now find that extra Si protrusions (black spheres) are induced in the \/7 x \/7 periodicity 
which are absent in the shorter periodicity. The protrusion generates the local tetrahedral 
geometry which is favorable for Si but at the same time induces stress energy around. 
The structure with the y/Y x ^/7 periodicity indicates that the incomplete protrusion is 
energetically unfavorable. Similarly, the AA- and slide-\/7 x \/7 structures also show the 
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extra protrusions. Such extra protrusions are also found in the 74^4-2x2 structure where 
four of 16 Si atoms are protruded in each cell while the AA-\pi x \/3 has one protruded 
atom in six Si atoms [Figs. [I](c) and [11(d)] . This excess protrusion is also observed in the 
highest-energy AB-lxl structure where half of Si atoms are protruded. The hex-OR-2x2 
structure which we have found has the lowest energy is free from such extra protrusion. 



FIG. 3. (Color online) Top view of OR-\/7 X \/7 structure. Red and green spheres are highly- 
protruded and moderately-protruded Si atoms as in Figs. [T^e) andl^f). Black spheres represent 
additionally protruded Si atoms. 

Figure 0] shows the PBE band structures of the four most stable bilayer silicene, slide- 
2x2, AA-\/Z X y/3, OR-\/3 x Cs and hex-OR-2x2. The total energy differences of these four 
structures are within the range of 31 meV/atom, indicating that they are stable at room 
temperature. The slide-2x2, AA-\f^ X \/3, and OR-\/3 x \/3 structures are narrow gap 
semiconductors with indirect band gaps of 0.07 (K-M line —>■ F), 0.06 (M ^ F), and 0.02 eV 
(M —7- M"), respectively. The band gaps of the three structures become approximately 0.37, 
0.34 and 0.40 eV, respectively, when we use the hybrid HSE functional. Interestingly, the 
indirect band gap (T-K line ^ K-M line) of the most stable hex-OR-2x2 structure is 0.8 eV 
(about 1.3 eV in the HSE calculation) considerably wide compared to other structures. 

The hex-OR-2x2 and OR-\/3 x \/3 structures have common protrusion but remarkably 
different band gaps. The conduction-band bottom state is more dispersive in the \/3 x 
s/Z structure compared to the 2x2 case. In fact, the Kohn-Sham wavefunction at the 
conduction-band bottom state of the OR-\/3 x structure have amplitude at unprotruded 



FIG. 4. PBE band structures of dynamically stable (a) slide-2 x 2, (b) AA-\pi> x \/3, (c) hex-OR- 
2x2, and (d) OR-\/3 x \/3 structures. The horizontal dashed line indicates the valence band top. 
The path is chosen so that valence band top and conduction-band bottom can be seen in the plot. 

Si atoms whereas the conduction-bottom state is highly-localized at the green protruded Si 
atoms in the 2x2 case (Fig. [5]). The deviation of the conduction-band bottom state can 
be related to the short distance between the protruded atoms in the \fZ x \fZ case. The 
band gap reduction is caused by the dispersive conduction band induced by the interaction 
between the protruded atoms. 

(a) (b) 




FIG. 5. (Color online) Conduction band bottom Kohn-Sham wavefunctions of (a) hex-OR-2x2 
and (b) OR-\/3 x \fZ structures. The yellow isosurfaces represent 1/3 of the maximum amplitude. 

In summary, we have performed density-functional calculations to clarify atomic and elec¬ 
tronic structures of freestanding bilayer silicene. After the extensive structural exploration, 
we have obtained four 1x1, two \/3 x \/3, and four 2x2 dynamically stable structures free 
from imaginary phonon frequencies. We have found that the atomic protrusion is a princi¬ 
pal structural characteristics which stabilizes bilayer silicene and induces optimum lateral 
periodicity. We have identihed the most stable bilayer silicene as the hex-OR-2x2 structure 
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and revealed that this is a semiconductor with the energy gap of 1.3 eV using the hybrid 
functional. 
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